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' We present a novel method for detecting nonlinearities, due to quantum electrodynamics through 

' photon-photon scattering, in Maxwell's equation. The photon-photon scattering gives rise to self- 

^SJ . interaction terms, which are similar to the nonlinearities due to the polarisation in nonlinear optics. 

These self-interaction terms vanish in the limit of parallel propagating waves, but if instead of 
parallel propagating waves the modes generated in waveguides are used, there will be a non-zero 
■ total effect. Based on this idea, we calculate the nonlinear excitation of new modes and estimate 

' the strength of this effect. Furthermore, we suggest a principal experimental setup. 
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According to QED, the non-classical phenomenon of photon-photon scattering can take place due to the exchange 
of virtual electron-positron pairs. The observation of this scattering is a long-standing experimental challenge which 



I \ has yet to be met. Photon-photon scattering is a second order effect (in terms of the fine structure constant a = 

^ ' /2eohc k, 1/137), and it can in standard notation be formulated in terms of the Euler-Heisenberg Lagrangian 

^ : density |l|,| 

^! ^ = eo^ + e(4^^ + 7^^) , (1) 

O ■ where 

^\ ^ ^ 45to4c5 ' 

^ = ^{E'^ — c^B"^), = cE • B, and and are the QED corrections. Here e is the electron charge, c the 
velocity of light, h the Planck constant and me the electron mass. We note that ^ = = in the limit of parallel 
propagating waves. The latter terms in ([T]) represent the effects of vacuum polarisation and magnetisation, and the 
_ QED corrected Maxwell's vacuum equations take the classical form using 

D^enE + P. H=— B M, 
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' where P and M arc of third order in the field amplitudes E and B, and /io = l/c^£o- Furthermore, they contain 
terms ^ and such that P = M = in the limit of parallel propagating waves. It is therefore necessary to 
use other waves in order to obtain an effect from these QED corrections. Several attempts have been presented in 
^ , the literature over the years |^|-^, where Refs. &3| mainly focused on principal issues, whereas the experimental 
. ^ possibilities for detection have been discussed in^]^. Soljacic & Segev concluded that, using their mechanism [Q, 
the detection of the QED nonlinearities will be technologically viable within ten to fifteen years, provided that the 
(-H laser power increases steadily. In this work we suggest the use of waveguides as an experimental setup, something 
O ^, which, as far as the authors know, has not been discussed previously for this purpose. The idea of using a waveguide 
' is to achieve a resonant coupling between the parallel propagating waves of different frequencies. We calculate the 
generated electromagnetic field for a rectangular waveguide, using the TEqi- and TEio-modes as pump waves. In the 
proposal for the experimental setup, the waveguide is replaced by a cavity, in order to prevent the convective loss 
of energy, and thereby maximise the saturated amplitude of the excited mode. The saturation level due to a finite 
conductivity of the cavity walls is estimated. We find that the excited mode can be detected even for moderate levels 
of the pump mode amplitudes, i.e. for field strengths that can be supported by the cavity walls. 

In a medium with polarisation P and magnetisation M the general wave equations for E and B are 
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V^B = /io V X (V X M) + —(V X P) 
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Furthermore, the effective polarisation and magnetisation in vacuum due to photon-photon scattering induced by the 
exchange of virtual electron-positron pairs are given by (see, e.g., Ref. [Q) 

P = 2^[2{E^ - c^B^)E + 7c^{E-B)B] , 

and 

M = -2c2^ [2{E^ - c^B^)B + 7(E • B)E] . 

Next we consider propagation in a rectangular waveguide with dimensions xq and yo (i-e- the region < a; < xq, 
< y < Uo is vacuum surrounded by walls that, as a starting point, are assumed to be perfectly conducting.). We 
assume that the TEiq- and TEgi-modes act as pump waves with distinct frequencies. To lowest order (i.e. neglecting 
the vacuum nonlinearities) the fields are 

Biz = Biz cos — exp[i(/ciz — uit)] + c.c. , (4a) 
Bix = - Biz sin ( — ] exp[i(fciz - ujit)] + c.c. , (4b) 
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Ely = Biz sm — exp[i(fciz - LOit)\ + c.c. , (4c) 



V TT / \Xo 

together with uj^ = kfc^ + tt^c^/xq for the TEio-mode, and 
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B2z = B2z cos — exp[i(fc2^ - ^^2^)] + c.c. , (5a) 

\yoJ 

-^2y = - (^~~^ sin (^-^^ exp[i(fc22; - ^2^)] + c.c. , (5b) 
E2x = - { ) B2z sin ( — ) exp[i(fc22: - ^2^)] + c.c. , (5c) 

V 7^ / \yaj 

together with = fc|c^ + it^c? jy^ for the TEoi-mode, where c.c. stands for complex conjugate. Here we have 
denoted the wave amplitudes, which to lowest order are constants, by Bzi and -Bz2, respectively. Substituting the 
linear expression for the fields into the cubic nonlinear terms, we note that there will be perturbations with frequency 
and wavenumber (cj3,/c3), where the possible combinations are (cj3,fc3) — (ioi^ki) , (0^2, fe) , (2loi ± uj2,2ki ± 
^2) , and (2ijj2 ± uji, 2^2 ± fci). If there is a small perturbation (e.g. of the order of 10~^^) of the amplitude of any of 
the original TE- modes, it would probably be a too difficult task to measure such an effect, whereas the appearance 
of a distinctly new frequency, although with small amplitude, will be easier to detect. We therefore concentrate on 
the two latter combinations, which are physically equivalent. Furthermore, if any of these combinations satisfies the 
dispersion relation for a natural mode of the waveguide, the amplitude of this resonantly driven mode will be much 
larger than the others. We therefore, for definiteness, choose to consider the matching condition 

([J3, fcs) = (2a;i - 072, 2fci - A;2) , (6) 

where the dimensions of the waveguide are assumed to be chosen such as to make (uj^jk^) a natural mode of the 
waveguide. 

Using the pump modes (^) and (^), we can get source terms in the wave equations either for a TEoi-mode or for 
a TMoi-mode. For the latter case, however, the frequency matching conditions and dispersions relations cannot be 
fulfilled simultaneously for real values of all wavenumbers, and thus we consider the excitation of a TEgi-mode, i.e. 
we let 
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Thus, using Eqs. (|^) and (|^), we can evaluate the source terms in Eqs. (|^) and (^). For a waveguide of finite length 
fulfilling the boundary condition of no incoming wave with frequency UJ3 at the waveguide starting point z = 0, we 
then have spatial growth. The ansatz B^z = B^z{z) exp[i(/c3Z — Lo^t)] -f c.c. in (^ thus gives a linear spatial growth 

B,z{z)^^^Bl^B*^, (7) 
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provided ^3 is not too small. Here the coupling constant V is 



£0 \ ^ ^J/o 



and the star in denotes complex conjugation. 

When designing parameters for an experiment, it might be tempting to choose parameters such that the dispersion 
relation is fulfilled for fcs = in which case we instead get a quadratic spatial growth of the amplitude. However, the 
case A:3 = is not the most interesting choice for two reasons. Firstly, the group velocity of and thereby the energy 
flux of the excited mode is proportional to fc3. Secondly it is impossible for a;3 > iO\^L02 to hold when ^3 = 0. The 
reason for requiring the excited wave to have a higher frequency than the others is that we then have the possibility to 
use waveguide filtering (see below) in order to measure the excited wave without any disturbing signals from the pump 
waves. Since k-j, cannot be small as compared to I/xq for > wi, a;2 to hold, the spatial growth in may be too slow 
for practical purposes, and we are therefore motivated to consider a cavity rather than a waveguide. Since the waves 
propagating in positive and negative directions are identical in that cavity, the boundary conditions imply temporal 
rather than spatial growth. If we assume that all waves have fc > (are propagating in the positive z-direction) 
in the waveguide example, the coupling coefhcient in a cavity can be found from the waveguide result, simply by 
noting that the positive propagating part of the standing pump waves couples to the positive propagating part of 
the excited standing wave, and vice versa. Since is a quadratic function of the wavenumbers, the same coupling 
strengths apply for standing waves ~ sin(n7rz/L) exp[— iwt], where L is the length of the cavity, as for propagating 
waves ^ exp[i(fcz — wt)], if we just let the wavenumber be fc = -izm: / L where the sign corresponds to the direction 
of propagation, and n is a positive integer. For design purposes we must keep in mind, however, that we now have 
additional^constraints relating the frequencies and dimensions since the wavenumbers cannot be chosen continously. 
If we let Bz represent the standing wave amplitude, and modify the ansatz such that the excited ampitude depends 
on time, we find that the temporal growth in a cavity is 

Bszit) = -z Bl^B^z . 

2^3 

Saturation occurs when the amplitude is large enough for linear damping due to a finite conductivity to balance the 
driving term. The saturated equilibrium amplitude is found to be 

|S3.cq| = |^|^L||B2,| , (8) 

where F is the linear damping rate of mode 3. If the linear damping is due to a finite conductivity cr of the walls, we 
have 

F = C/i^ 

in the regime F/W3 ^ 1, where 5 = Re (i/^oO'W3)^^^ is the skindepth, R = Re((T^^)/(5 is the surface resistance of the 

cavity walls, and C/ is a function of the geometry that is of the order of unity. 

In order to avoid problems due to the large amplitude pump signals when trying to measure the excited mode, 
it is convenient to use a slightly modified cavity rather than an idealized cubic cavity. Such a cavity could consist 
of two parts; one cavity (I), in which the signal is generated, attached to another cavity (II), acting as a waveguide 
filter. The dimensions of cavity I should be chosen such as to keep the frequency of the excited mode above cut-off in 
cavity II, whereas the pump modes should be below cut-off. By measuring the excited signal far enough into cavity 
II (in practice 30-40 pump decay lengths), the pump waves effectively vanish, and we can forget about the disturbing 
influence of the pump waves. However, note that a certain fine tuning of the length of cavity II might be necessary 
to keep a maximally efficient phase relation between the excited eigenmode and the pump waves in cavity I. 

Next we demonstrate that the experimental model setup gives signals that can be detected with presently available 
technology. High performance, i.e. large electromagnetic fields combined with low dissipative losses, can be found 
in superconducting cavities, which among other things are used for particle accelerator purposes Adopting data 
from these experiments, we assume that the pump waves have a field-strength Ecnt ~ 30MV/m, i.e. close to the 
maximum that can be tolerated by the walls without field emissions. For a cavity resistance R ^ \ nfi, corresponding 
to superconducting niobium at a temperature 1.4 K and a frequency u;3 ~ 2 x lO"'^" rad/s, we find from Eq. (||) that 
the saturated energy flux P3 of mode 3 is of the order of P3 ~ 10~^ W/m^. Here we have used the simple order of 
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magnitude estimate that all wavelengths are comparable to the dimensions Xq and yo when evaluating V. Clearly 
this energy flux is above the detection level by several orders of magnitude. Note, however, the importance of the 
superconducting walls for the output level of the excited mode. For copper at room temperature, the cavity resistance 
increases by a factor ~ 10^ as compared to the above example, and consequently the energy flux of the excited mode 
falls by a factor ~ 10"^''. In this latter case it is questionable whether the excited signal can be detected. 

To our knowledge, it has not been possible to verify the Euler-Heisenberg Lagrangian density experimentally. 
The above calculation of the QED mode coupling strength and the subsequent estimations suggest that it can be 
very suitable to use two pump modes in a superconducting cavity for this purpose. The parameters of the problem 
should be designed such as to simultaneously fulfill the dispersion relations of each mode together with the matching 
condition (^. Naturally, care must be taken when drawing the conclusions, since there are certain effects that we have 
not yet addressed. For example, in the present model the conductivity of the walls is linear, but in principle there 
might be a nonlinear contribution to the conductivity that could give rise to a small signal at the same frequency 
as the QED contribution. On the other hand, to our knowledge, there are no theoretical or experimental reports of 
such effects. Secondly, in reality the vacuum in the waveguide will not be perfect, and in principle this may lead to 
dielectric breakdown for the pump field strength considered in the above estimate. However, we do not expect this 
particular effect to be a serious threat against our proposal, since similar electric field strengths have been reached 
in present experiments 1^, and the pump field strength that is actually needed is much less than the one considered 
in the estimate. Thus we conclude that it is likely that the effect of photon-photon scattering in vacuum due to the 
exchange of virtual electron-positron pairs can be meaured using existing technology. 
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